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A programmed computer is used to create 
a digital model of an individual component of a 
patients dentition. The computer obtains a 3D 
digital model of the patient's dentition, identities 
points (766) in the dentition model that lie on an 
inter-proximal margin between adjacent teeth in the 
patient's dentition, and uses the identified points to 
create a cutting surface (782) that separates portions 
of the dentition model representing the adjacent 
teeth. 
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MANIPULATING A DIGITAL DENTITION MODEL TO FORM MODELS OF 
INDIVIDUAL DENTITION COMPONENTS 



10 



15 

TECHNOLOGICAL FIELD 

The invention relates to the fields of computer-assisted dentistry and orthodontics. 

20 BACKGROUND 

Two-dimensional (2D) and three-dimensional (3D) digital image technology has recently 
been tapped as a tool to assist in dental and orthodontic treatment. Many treatment providers use 
some form of digital image technology to study the dentitions of patients. U.S. patent 
application 09/169,276 describes the use of 2D and 3D image data in forming a digital model of 
25 a patient's dentition, including models of individual dentition components. Such models are 

useful, among other things, in developing an orthodontic treatment plan for the patient, as well as 
in creating one or more orthodontic appliances to implement the treatment plan. 

SUMMARY 

30 The inventors have developed several computer-automated techniques for subdividing, or 

segmenting, a digital dentition model into models of individual dentition components. These 
dentition components include, but are not limited to, tooth crowns, tooth roots, and gingival 
regions. The segmentation techniques include both human-assisted and fully-automated 
techniques. Some of the human-assisted techniques allow a human user to provide "algorithmic 

35 hints" by identifying certain features in the digital dentition model. The identified features then 
serve as a basis for automated segmentation. Some techniques act on a volumetric 3D image 
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model, or "voxel representation," of the dentition, and other techniques act on a geometric 3D 
model, or "geometric representation." 

In one aspect, a computer implementing the invention receives a data set that forms a 
three-dimensional (3D) representation of the patient's dentition, applies a test to the data set to 

5 identify data elements that represent portions of the individual component, and creates a digital 
model of the individual component based upon the identified data elements. Some 
implementations require the computer to identify data elements that form one or more 2D cross- 
sections of the dentition in one or more 2D planes intersecting the dentition. In many of these 
embodiments, these 2D planes are roughly parallel to the dentition's occlusal plane. The 

10 computer analyzes the features of the 2D cross-sections to identify data elements that correspond 
to the individual component to be modeled. For example, one technique requires the computer 
to identify cusps in the 2D cross-sectional surface of the dentition, where the cusps represent the 
locations of an interproximal margin between teeth in the dentition. One variation of this 
technique allows the computer to confine its search for cusps in one 2D plane to areas in the 

1 5 vicinity of cusps already identified on another 2D plane. Another variation allows the computer 
to link cusps on adjacent 2D planes to form a solid surface representing the interproximal 
margin. Some embodiments allow the computer to receive input from a human user identifying 
the cusp locations in one or more of the 2D cross sections. 

Other embodiments require the computer to identify data elements that represent a 

20 structural core, or skeleton, of each individual component to be modeled. The computer creates 
the model by linking other data elements representing the individual component to the structural 
core. 

In another aspect, a computer implementing the invention receives a three-dimensional 
(3D) data set representing the patient's dentition, applies a test to identify data elements tHat 

25 represent an interproximal margin between two teeth in the dentition, and applies another 
computer-implemented test to select data elements that lie on one side of the interproximal 
margin for inclusion in the digital model. Some implementations require the computer to 
identify data elements that form one or more 2D cross-sections of the dentition in one or more 
2D planes intersecting the dentition roughly parallel to the dentition's occlusal plane. 

30 In another aspect, a computer implementing the invention receives a 3D data set 

representing at least a portion of the patient's dentition, including at least a portion of a tooth and 
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gum tissue surrounding the tooth; applies a test to identify data elements lying on a gingival 
boundary that occurs where the tooth and the gum tissue meet; and applies a test to the data 
elements lying on the boundary to identify other data elements representing portions of the tooth. 
In one aspect, the invention involves obtaining a three-dimensional (3D) digital model of 

5 a patient's dentition and analyzing the model to determine the orientation of at least one axis of 
the model automatically. In some implementations, the model's z-axis is found by creating an . 
Oriented Bounding Box (OBB) around the model and identifying the direction in which the OBB 
has minimum thickness. The z-axis extends in this direction, from the model's bottom surface to 
its top surface. Moreover, in a dentition model having only one mandible, one of the model 

10 surfaces is substantially flat and an opposite surface is textured. The direction of the positive z- 
axis can be identified in this type of model by identifying which of the surfaces is flat or 
textured. One technique for doing so involves creating one or more planes that are roughly 
normal to the z-axis and then creating line segments that extend between the planes and the top 
and bottom surfaces of the dentition model. The surface for which all of the line segments are of 

15 one length is identified as being the flat surface, and the surface for which the line segments have 
varying lengths is identified as being the textured surface. 

In other implementations, the x- and y-axes are found by selecting a two-dimensional 
(2D) plane that contains the axes and an arch-shaped cross section of the dentition model and 
identifying the orientations of the axes in this plane. In general, the arch-shaped cross section is 

20 roughly symmetrical about the y-axis. One technique for identifying the y-axis involves 
identifying a point at each end of the arch-shaped cross section, creating a line segment that 
extends between the identified points, and identifying the orientation of the y-axis as being 
roughly perpendicular to the line segment. The point at each end of the arch can be identified by 
selecting a point that lies within an area surrounded by the arch-shaped cross section, creating a 

25 line segment that extends between the selected point and an edge of the 2D plane, sweeping the 
line segment in a circular manner around the selected point, and identifying points at the ends of 
the arch-shaped cross section at which the sweeping line segment begins intersecting the cross 
section of the dentition model and stops intersecting the cross section of the dentition model. In 
general, the x-axis is perpendicular to the y-axis. 

30 In another aspect, the invention involves using a programmed computer to create a digital 

model of an individual component of a patient's dentition by obtaining a 3D digital model of the 
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patient's dentition, identifying points in the dentition model that lie on an inter-proximal margin 
between adjacent teeth in the patient's dentition, and using the identified points to create a cutting 
surface for use in separating portions of the dentition model representing the adjacent teeth. 
In some implementations, 2D cross sections of the dentition model are displayed to a 
5 human operator, and the operator provides input identifying approximate points at which the 
interproximal margin between the adjacent teeth meets gingival tissue. In some cases, the 
dentition model includes a 3D volumetric model of the dentition, and the input provided by the 
operator identifies two voxels in the volumetric model. The computer then defines a 
neighborhood of voxels around each of the two voxels identified by the human operator, where 

1 0 each neighborhood includes voxels representing the dentition model and voxels representing a 
background image. The computer selects the pair of voxels, one in each neighborhood, 
representing the background image that lie closest together. 

In some of these implementations, the computer also identifies voxels on another 2D 
cross section that represent the interproximal margin. One technique for doing so is by defining 

15 a neighborhood of voxels around each of the selected voxels, where each neighborhood includes 
voxels representing the dentition model and voxels representing a background image, projecting 
the neighborhoods onto the other 2D cross section, and selecting two voxels in the projected 
neighborhoods that represent the inter-proximal margin. 

In another aspect, the invention involves displaying an image of a dentition model, 

20 receiving input from a human operator identifying points in the image representing a gingival 
line at which a tooth in the dentition model meets gingival tissue, and using the identified points 
to create a cutting surface for use in separating the tooth from the gingival tissue in the dentition 
model. The cutting surface often extends roughly perpendicular to the dentition's occlusal plane. 
In some implementations, the cutting surface is created by projecting at least a portion of 

25 the gingival line onto a plane that is roughly parallel to the occlusal plane and then creating a 
surface that connects the gingival line to the projection. One way of establishing the plane is by 
fitting the plane among the points on the gingival line and then shifting the plane away from the 
tooth in a direction that is roughly normal to the plane. For example, the plane can be shifted 
along a line segment that includes a point near the center of the tooth and that is roughly 

30 perpendicular to the plane. The length of the line segment usually approximates the length of a 
tooth root. 
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In other embodiments, the cutting surface extends roughly parallel to the dentition's 
occlusal plane in the dentition model. In some of these embodiments, the input received from 
the human operator identifies points that form two 3D curves representing gingival lines at which 
teeth in the dentition model meet gum tissue on both the buccal and lingual sides of the dentition 

5 model. The cutting surface is created by fitting a surface among the points lying on the two 
curves. For each tooth, a point lying between the two curves is identified and surface triangles . 
are created between the identified point and points on the two curves. One technique for 
identifying the point involves averaging, for each tooth, x, y and z coordinate values of the points 
on portions of the two curves adjacent to the tooth, 

1 0 Other embodiments involve creating, for each tooth, a surface that represents the tooth's 

roots. One technique for doing so involves projecting points onto a plane that is roughly parallel 
to the occlusal plane and connecting points on the two curves to the projected points. The 
surface can be used to separate portions of the dentition model representing the tooth roots from 
portions representing gingival tissue. The model of the tooth roots is then connected to die tooth 

15 model. 

Other embodiments and advantages are apparent from the detailed description and the 
claims below. 

DESCRIPTION OF THE DRAWINGS 

20 FIGS. 1A, IB, and 2 are partial views of a dentition model as displayed on a computer 

monitor and segmented with a human-operated saw tool. 

FIG. 3 is a partial view of a dentition model as displayed on a computer monitor and 
segmented with a human-operated eraser tool. 

FIG. 4 is a view of a dentition model for which a feature skeleton has been identified. 
25 FIGS. 5 and 6 are flowcharts for a feature skeleton analysis technique used in segmenting 

a dentition model. 

FIG. 7A is a horizontal 2D cross-sectional view of a dentition model 

FIG. 7B is a side view of a dentition model intersected by several 2D planes. 

FIG. 8 is a flowchart for a 2D slice analysis technique used in segmenting a dentition 

30 model. 
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FIGS. 9 and 10A through 10C each shows a group of voxels in a 2D slice of a dentition 

model. 

FIG. 1 1 is a flowchart for an automatic cusp detection technique used in segmenting a 
dentition model. 

FIG. 12 is a horizontal 2D cross section of a dentition model illustrating a neighborhood 
filtered automatic cusp detection technique used in segmenting the dentition model. 

FIG. 13 shows two groups of voxels in a 2D slice of a dentition model illustrating the 
neighborhood filtered automatic cusp detection technique. 

FIG. 14 is a flowchart for the neighborhood filtered automatic cusp detection technique. 

FIG. 15 is a horizontal 2D cross section of a dentition model illustrating an arch curve 
fitting technique used in segmenting the dentition model. 

FIG. 16 is a flowchart for the arch curve fitting technique. 

FIG. 17 is a horizontal 2D cross section of a dentition model illustrating a curve creation 
technique for use with the arch curve fitting technique. 

FIG. 18 is a flowchart for the curve creation technique. 

FIGS. 19A and 19B are a perspective view and a vertical 2D cross-sectional view of a 
dentition model illustrating another technique for use in segmenting the dentition model. 

FIGS. 20 and 21 are flowcharts of the technique illustrated in FIGS. 19A and 19B. 

FIG. 22 is a vertical 2D cross-sectional view of a dentition model illustrating the gingival 
margin detection technique for use in segmenting the dentition model. 

FIG. 23 shows a group of voxels in a 2D slice of a dentition model illustrating a gingival 
margin detection technique. 

FIG. 24 is a flowchart for the gingival margin detection technique. 

FIG. 25 shows a digital dentition model inside an Oriented Bounding Box (OBB). 

FIG. 26 illustrates a technique for properly orienting a digital dentition model along a 

z-axis. 

FIGS. 27 A, 27B, and 27C illustrate a technique for properly orienting a digital dentition 
model along x- and y-axes. 

FIGS. 28, 29, 30 and 31 are flowcharts for the techniques of FIGS. 25, 26, and 27A-C. 
FIGS. 32 and 33 illustrate a human-assisted technique for identifying interproximal 
margins between teeth. 
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FIG. 34 is a flowchart for the technique of FIGS. 32 and 33. 

FIGS. 35A through 35F illustrate a technique for segmenting a digital dentition model 
into models of individual teeth and gum tissue. 

FIG. 36 is a flowchart for the technique of FIGS. 35A through 35F. 
5 FIGS. 37A, 37B, and 37C illustrate another technique for segmenting a digital dentition 

model into models of individual teeth. 

FIGS. 38 and 39 are flowcharts for the technique of FIGS. 37A, 37B, and 37C. 

DETAILED DESCRIPTION 

1 0 U.S. patent application 09/1 69,276 describes techniques for generating a 3D digital data 

set that contains a model of a patient's dentition, including the crowns and roots of the patient's 
teeth as well as the surrounding gum tissue. One such technique involves creating a physical 
model of the dentition from a material such as plaster and then digitally imaging the model with 
a laser scanner or a destructive scanning system. These techniques are used to produce a digital 

15 volumetric 3D model ("volume element representation" or "voxel representation") of the 
dentition model, and/or a digital geometric 3D surface model ("geometric model") of the 
dentition. The computer-implemented techniques described below act on one or both of these 
types of 3D dentition models. 

In creating a voxel representation, the physical model is usually embedded in a potting 

20 material that contrasts sharply with the color of the physical model to enhance detection of the 
dentition features. A white dentition model embedded in a black potting material provides the 
sharpest contrast A wide variety of information can be used to enhance the 3D model, including 
data taken from photographic images, 2D and 3D x-rays scans, computed tomography (CT) 
scans, and magnetic resonance imaging (MRI) scans of the patient's dentition. 

25 The 3D data set is loaded into a computer which, under control of a program 

implementing one or more techniques of the dentition, either with or without human assistance, 
segments the digital dentition model into digital models of individual dentition components, such 
as teeth and gingival tissue. In one implementation, the computer produces a digital model of 
each individual tooth in the patient's dentition, as well as a digital model of the gingival tissue 

30 j surrounding the teeth. 
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To segment the digital dentition model accurately, the computer often must know the 
exact orientation of the dentition model. One technique for establishing the orientation of the 
digital dentition model in the 3D data set involves holding the physical dentition model at a 
prescribed orientation during the digital imaging process discussed above. Embedding the 
5 physical model at a particular orientation in a solid potting material is one way of holding the 
physical model. In some systems, however, even this technique introduces small errors in the 
orientation of the dentition model. 

Orienting the- Digital Dentition Model 

10 FIGS. 25, 26, 27A-C and 28 illustrate several techniques used by the computer to orient 

the digital dentition model 500 properly. The computer first obtains a digital model of the 
dentition using one of the techniques described above (step 700). The computer then locates the 
model's z-axis 502, which in the depicted example extends from the base of the model toward the 
roof of the patient's mouth and is normal to the dentition's occlusal plane (step 702). The 

15 computer then locates the models y-axis 504, which in the depicted example extends from an 
area lying within the dental arch toward the patient's front teeth (step 704). Using the right-hand 
rule, the computer then defines the model's x-axis 506 to extend from an area lying within the 
dental arch toward the teeth on the right side of the patient's mouth (step 706). The occlusal 
plane is a plane that is pierced by all of the cusps of the patient's teeth when the patient's 

20 mandibles interdigitate. Techniques for identifying the occlusal plane include receiving user 
input identifying the location of the plane and conducting a fully-automated analysis of the 
dentition model. 

FIGS. 25, 26, and 29 show one technique for identifying the z-axis 502. The computer 
first identifies the dentition model 500 in the 3D data set (step 710). For 3D geometric data, 

25 identifying the dentition model is simply a matter of locating the geometric surfaces. For 3D 

volumetric data, identifying the dentition model involves distinguishing the lighter voxels, which 
represent the dentition model, from the darker voxels, which represent the background. The 
computer then fits an Oriented Bounding Box ("OBB") 510 around the dentition model 500 
using a conventional OBB fitting technique (step 712). The dimension in which the OBB 510 

30 has its smallest thickness Tmin is the dimension in which the z-axis 502 extends (step 714). 
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After determining the dimension in which the z-axis extends 502, the computer 
determines whether the dentition model is facing upward or downward, z.e., in which direction 
the positive z-axis extends. FIGS. 26 and 30 illustrate a technique for determining the direction 
of the positive z-axis. This technique relies on an observation that the bottom surface 512 of the 
5 dentition model is flat and the upper surface 514 follows the jagged contours of the patient's 
teeth. This technique also relies on an assumption that the model at this point includes only one 
of the patient's mandibles. 

The computer first creates one or more planes 516, 518 that are normal to the z-axis 502 
(step 720). The computer then creates line segments 515A, 515B between the planes 516, 518 

1 0 and the surfaces 5 12, 5 14 of the model (step 722). The line segments 5 1 5A that touch the flat 
bottom surface 512 are all of approximately the same length (step 724). The line segments 515B 
that touch the jagged top surface 514 have varying lengths (step 726). The computer identifies 
the positive z-axis as extending from the bottom surface 512 toward the top surface 514 and 
orients the digital dentition model 500 accordingly (step 728). 

15 FIGS. 27A, 27B, 27C, and 31 illustrate a technique for identifying the y-axis 504 and the 

x-axis 506 of the dentition model 500. The computer begins by selecting a 2D slice 520 of data 
that is normal to the z-axis and that contains a cross section 522 of the dentition model 
(step 730). This technique relies on an observation that the cross section 522 of the dentition 
model is arch shaped. The computer identifies a point 524 at or near the center of the 2D slice 

20 520 (step 732). The computer then creates a line segment 526 (or 530) that extends from the 

selected point 524 to an edge 528 (or 532) of the slice 520 (step 734). The direction in which the 
line segment extends is arbitrary, so the line segment may or may not intersect the dental cross 
section. The depicted example shows two line segments 526, 530, one of which intersects the 
dental cross section 522, the other of which does not. 

25 The computer then begins rotating, or sweeping, one of the line segments 526, 530 about 

the center point 524 (step 736). In general, the computer sweeps the line segment in small, 
discrete steps, usually on the order of five degrees of rotation. As it is swept, a line segment 526 
that initially intersects the dental cross section 522 will eventually stop intersecting the cross 
section 522, and the computer marks the point 534 at which this occurs. As sweeping continues, 

30 the line segment 526 will eventually resume intersecting the cross section 522, and the computer 
marks the point 536 at which this occurs. Likewise, a line segment 530 that initially does not 
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intersect the cross section 522 eventually will begin intersecting the cross section 522, and the 
computer marks the point 536 at which this occurs. The computer also marks the point 534 at 
which this line segment 530 stops intersecting the cross section 522 (step 738). The computer 
stops sweeping the line segments 526, 530 after marking both of the points 534, 536 (step 740). 
5 The computer then creates a line segment 538 that extends between the two marked 

points 534, 536 (step 742). The y-axis 504 of the dentition model extends roughly normal to this 
line segment 538 through the front 540 of the dental arch (step 744). The x-axis 506 extends 
roughly parallel to this line segment 538 through the right side 542 of the dental arch (step 746). 
The computer uses this line segment 538 to orient the dentition model correctly along the x- and 
10 y-axes (step 748). 

Segmenting the Digital Dentition Model Into Individual Component Models 
Some computer-implemented techniques for segmenting a 3D dentition model into 
models of individual dentition components require a substantial amount of human interaction 

15 with the computer. One such technique, which is shown in FIGS. 1 A, IB, and 2, provides a 
graphical user interface with a feature that imitates a conventional saw, allowing the user to 
identify components to be cut away from the dentition model 100. The graphical user interface 
provides a rendered 3D image 100 of the dentition model, either at one or more static views from 
predetermined positions, as shown in FIGS. 1 A and IB, or in a "full 3D" mode that allows the 

20 user to alter the viewing angle, as shown in FIG. 2. The saw tool is implemented as a set of 
mathematical control points 102, represented graphically on the rendered image 100, which 
define a 3D cutting surface 104 that intersects the volumetric or geometric dentition model. The 
computer subdivides the data elements in the dentition model by performing a surface 
intersection operation between the 3D cutting surface 104 and the dentition model. The user sets 

25 the locations of the mathematical control points, and thus the geometry and position of the 3D 
cutting surface, by manipulating the control points in the graphical display with an input device, 
such as a mouse. The computer provides a visual representation 104 of the cu_ing surface on 
the display to assist the user in fitting the surface around the individual component to be 
separated. Once the intersection operation is complete, the computer creates a model of the 

30 individual component using the newly segmented data elements. 
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Another technique requiring substantial human interaction, shown in FIG. 3, is a 
graphical user interface with a tool that imitates a conventional eraser. The eraser tool allows the 
user to isolate an individual dentition component by removing portions of the dentition model 
that surround the individual component. The eraser tool is implemented as a 3D solid 110, 
5 typically having the shape of a rectangular prism, or a curved surface that matches the shape of a 
side surface of a tooth. The solid is made as small as possible, usually only a single voxel thick, 
to minimize degradation of the data set As with the saw technique above, the graphical user 
interface presents the user with a rendered 3D image 1 12 of the dentition model at one or more 
predetermined static views or in a full 3D mode. The user identifies portions of the dentition 

1 0 model for removal by manipulating a graphical representation 1 1 0 of the 3D solid with an input 
device. In alternative embodiments, the computer either removes the identified portions of the 
dentition model as the user moves the eraser 1 12, or the computer waits until the user stops 
moving the eraser and provides an instruction to remove the identified portions. The computer 
updates the display in real time to show the path 1 14 of the eraser through the dentition model. 

15 Other computer-implemented segmentation techniques require little or no human 

interaction during the segmentation process. One such technique, which is illustrated in FIG. 4, 
involves the application of conventional "feature skeleton" analysis to a volumetric 
representation of the dentition model. This technique is particularly useful in identifying and 
modeling individual teeth. In general, a computer applying this technique identifies a core of 

20 voxels, that forms a skeleton 122 for the dentition 120. The skeleton 122 roughly resembles the 
network of biological nerves within patient's teeth. The computer then divides the skeleton 122 
into branches 124, each containing voxels that lie entirely within one tooth. One technique for 
identifying the branches is by defining a plane 126 that cuts through the skeleton 122 roughly 
parallel to the occlusal plane of the patient's dentition (horizontal plane"). Each branch 124 

25 intersects the horizontal plane 126 at one or more points, or clusters, that are relatively distant 
from the clusters associated with the other branches. The computer forms the individual tooth 
models by linking other voxels to the appropriate branches 124 of the skeleton. 

FIG. 5 describes a particular technique for forming a skeleton in the dentition model. 
The computer first identifies the voxels in the dentition model that represent the tooth surfaces 

30 (step 130). For a voxel representation that is created from a physical model embedded in a 

sharply contrasting material, identifying the tooth surfaces is as simple as identifying the voxels 
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at which sharp changes in image value occur, as described in U.S. patent application 09/169,276. 
The computer then calculates, for each voxel in the model, a distance measure indicating the 
physical distance between the voxel and the nearest tooth surface (step 132). The computer 
identifies the voxels with the largest distance measures and labels each of these voxels as 
5 forming a portion of the skeleton (step 134). Feature skeleton analysis techniques are described 
in more detail in the following publications: (1) Gagvani and Silver, "Parameter Controlled 
Skeletons for 3D Visualization," Proceedings of the IEEE Visualization Conference (1997); 
(2) Bertrand, "A Parallel Thinning Algorithm for Medial Surfaces," Pattern Recognition Letters, 
v. 16, pp. 979-986 (1995); (3) Mukherjee, Chatterji, and Das, "Thinning of 3-D Images Using 
10 the Safe Point Thinning Algorithm (SPTA)," Pattern Recognition Letters, v. 10, pp. 167-173 
(1989); (4) Niblack, Gibbons, and Capson, "Generating Skeletons and Centerlines from the 
Distance Transform," CVGIP: Graphical Models and Image Processing, v. 54, n. 5, pp. 420-437 
(1992). 

Once a skeleton has been formed, the computer uses the skeleton to divide the dentition 

15 model into 3D models of the individual teeth. FIG. 6 shows one technique for doing so. The 
computer first identifies those portions of the skeleton that are associated with each individual 
tooth. To do so, the computer defines a plane that is roughly parallel to the dentition's occlusal 
surface and that intersects the skeleton near its base (step 136). The computer then identifies 
points at which the plane and the skeleton intersect by identifying each voxel that lies on both the 

20 skeleton and the plane (step 138). In general, a single tooth includes all of the voxels that lie in a 
particular branch of the skeleton; and because the plane intersects the skeleton near its base, 
voxels that lie together in a branch of the skeleton usually cluster together on the intersecting 
plane. The computer is able to locate the branches by identifying voxels on the skeleton that lie 
within a particular distance of each other on the intersecting plane (step 140). The computer then 

25 identifies and labels all voxels on the skeleton that belong to each branch (step 142). 

Once the branches are identified, the computer links other voxels in the model to the 
branches. The computer begins by identifying a reference voxel in each branch of the skeleton 
(step 144). For each reference voxel, the computer selects an adjacent voxel that does not lie on 
the skeleton (step 146). The computer then processes the selected voxel, determining whether 

30 the voxel lies outside of the dentition, z.e., whether the associated image value is above or below 
a particular threshold value (step 148); determining whether the voxel already is labeled as 
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belonging to another tooth (step 150); and determining whether the voxel's distance measure is 
greater than the distance measure of the reference voxel (step 1 52). If none of these conditions is 
true, the computer labels the selected voxel as belonging to the same tooth as the reference voxel 
(step 1 54). The computer then repeats this test for all other voxels adjacent to the reference 
5 voxel (step 156). Upon testing all adjacent voxels, the computer selects one of the adjacent 
voxels as a new reference point, provided that the adjacent voxel is labeled as belonging to the 
same tooth, and then repeats the test above for each untested voxel that is adjacent to the new 
reference point This process continues until all voxels in the dentition have been tested. 

FIGS. 7A and 7B illustrate another technique for identifying and segmenting individual 

10 teeth in the dentition model. This technique, called "2D slice analysis," involves dividing the ' 
voxel representation of the dentition model into a series of parallel 2D planes 160, or slices, that 
are each one voxel thick and that are roughly parallel to the dentition's occlusal plane, which is 
roughly normal to the model's z-axis. Each of the 2D slices 160 includes a 2D cross section 162 
of the dentition, the surface 1 64 of which represents the lingual and buccal surfaces of the 

15 patient's teeth and/or gums. The computer inspects the cross section 162 in each 2D slice 160 to 
identify voxels that approximate the locations of the interproximal margins 166 between the 
teeth. These voxels he at the tips of cusps 165 in the 2D cross-sectional surface 164. The 
computer then uses the identified voxels to create 3D surfaces 168 intersecting the dentition 
model at these locations. The computer segments the dentition model along these intersecting 

20 surfaces 168 to create individual tooth models. 

FIG. 8 describes a particular implementation of the 2D slice analysis technique. The 
computer begins by identifying the voxels that form each of the 2D slices (step 170). The 
computer then identifies, for each 2D slice, the voxels that represent the buccal and lingual 
surfaces of the patient's teeth and gums (step 172) and defines a curve that includes all of these 

25 voxels (step 174). This curve represents the surface 164 of the 2D cross section 162. 

The computer then calculates the rate of curvature (i.e., the derivative of the radius of 
curvature) at each voxel on the 2D cross-sectional surface 164 (step 176) and identifies all of the 
voxels at which local maxima in the rate of curvature occur (step 178). Each voxel at which a 
local maximum occurs represents a "cusp" in the 2D cross-sectional surface 164 and roughly 

30 coincides with an interproximal margin between teeth. In each 2D slice, the computer identifies 
pairs of these cusp voxels that correspond to the same interproximal margin (step 180), and the 
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computer labels each pair to identify the interproximal margin with which it is associated (step 
182). The computer then identifies the voxel pairs on all of the 2D slices that represent the same 
interproximal margins (step 184). For each interproximal margin, the computer fits a 3D surface 
168 approximating the geometry of the interproximal margin among the associated voxel pairs 
5 (step 186). 

FIG. 9 illustrates one technique for creating the 3D surfaces that approximate the 
interproximal margins. For each pair of cusp voxels 190a-b in a 2D slice that are associated with 
a particular interproximal region, the computer creates a line segment 192 bounded by these cusp 
voxels 190a-b. The computer changes the colors of the voxels in the line segment, including the 

10 cusp voxels 190a-b that bound the segment, to contrast with the other voxels in the 2D slice. The 
computer creates line segments in this manner in each successive 2D slice, forming 3D surfaces 
that represent the interproximal regions. All of the voxels that lie between adjacent ones of these 
3D surfaces represent an individual tooth. 

FIGS. 10A through 10C illustrate a refinement of the technique shown in FIG. 9. The 

15 refined technique involves the projection of a line segment 200 from one slice onto a line 

segment 206 on the next successive slice to form, for the associated interproximal margin, a 2D 
area bounded by the cusp voxels 202a-b, 204a-b of the line segments 200, 206. If the line 
segments 200, 206 are oriented such that any voxel on one segment 200 is not adjacent to a voxel 
on the other segment 206, as shown in FIG. 10A, then the resulting 3D surface is discontinuous, 

20 leaving unwanted "islands" of white voxels 208, 2 1 0. 

The computer eliminates these discontinuities by creating two new line segments 212, 
214, each of which is bounded by one cusp voxel 202a-b, 204a-b from each original line segment 
200, 206, as shown in FIG. 10B. The computer then eliminates the islands between the new line 
segments 212, 214 by changing the colors of all voxels between the new line segments 212, 214, 

25 as shown in FIG. 1 0C. 

Automated segmentation is enhanced through a technique known as "seed cusp 
detection." The term "seed cusp" refers to a location at which an interproximal margin between 
adjacent teeth meets the patient's gum tissue. In a volumetric representation of the patient's 
dentition, a seed cusp for a particular interproximal margin is found at the cusp voxel that lies 

30 closest to the gum line. By applying the seed cusp detection technique to the 2D slice analysis, 
the computer is able to identify all of the seed cusp voxels in the 3D model automatically. 
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FIG. 1 1 shows a particular implementation of the seed cusp detection technique, in which 
the computer detects the seed cusps by identifying each 2D slice in which the rate of curvature of 
a cusp first falls below a predetermined threshold value. The computer begins by selecting a 2D 
slice that intersects all of the teeth in the arch (step 220). The computer attempts to select a slice 

5 that is near the gingival regions but that does not include any voxels representing gingival tissue. 
The computer then identifies all of the cusp voxels in the 2D slice (step 222). If the rate of 
curvature of the 2D cross section at any of the cusp voxels is less than a predetermined threshold 
value, the computer labels that voxel as a seed cusp (step 224). The computer then selects the 
next 2D slice, which is one voxel layer closer to the gingival region (step 226), and identifies all 

10 of the cusp voxels that are not associated with a cusp for which the computer has already 

identified a seed cusp (step 228). If the rate of curvature of the 2D cross section is less than the 
predetermined threshold value at any of these cusp voxels, the computer labels the voxel as a 
seed cusp (step 230) and proceeds to the next 2D slice. The computer continues in this manner 
until a seed cusp voxel has been identified for each cusp associated with an interproximal margin 

15 (step 232). 

FIGS. 32, 33, and 34 illustrate a human-assisted technique, known as "neighborhood- 
filtered seed cusp detection," for detecting seed cusps in the digital dentition model. This 
technique allows a human operator to scroll through 2D image slices on a video display and 
identify the locations of the seed cusps for each of the interproximal margins. The computer 

20 displays the 2D slices (step 750), and the operator searches the 2D slices to determine, for each 
adjacent pair of teeth, which slice 550 most likely contains the seed cusps for the corresponding 
interproximal margin. Using an input device such as a mouse or an electronic pen, the user 
marks the points 552, 554 in the slice 550 that appear to represent the seed cusps (step 752). 
With this human guidance, the computer automatically identifies two voxels in the slice as the 

25 seed cusps. 

The points 552, 554 identified by the human operator may or may not be the actual seed 
cusps 560, 562, but these points 552, 554 lie very close to the actual seed cusps 560, 562. As a 
result, the computer confines its search for the actual seed cusps 560, 562 to the voxel 
neighborhoods 556, 558 immediately surrounding the points 552, 554 selected by the human 
30 operator. The computer defines each of the neighborhoods 556, 558 to contain a particular 
number of voxels, e.g., twenty-five arranged in a 5 x 5 square, as shown here (step 754). The 
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computer then tests the image values for all of the voxels in the neighborhoods 556, 558 to 
identify those associated with the background image and those associated with the dentition (step 
756). In this example, voxels in the background are black and voxels in the dentition are white. 
The computer identifies the actual seed cusps 560, 562 by locating the pair of black voxels, one 
from each of the neighborhoods 556, 558, that lie closest together (step 758). In the depicted 
example, each of the actual seed cusps 560, 562 lies next to one of the points 552, 554 selected 
by the human operator. 

FIGS. 12, 13, and 14 illustrate a technique, known as "neighborhood-filtered cusp 
detection," by which the computer focuses its search for cusps on one 2D slice to neighborhoods 
244, 246 of voxels defined by a pair of previously detected cusp voxels 240, 242 on another 2D 
slice. This technique is similar to the neighborhood-filtered seed cusp detection technique 
described above. 

Upon detecting a pair of cusp voxels 240, 242 in a 2D slice at level N (step 250), the 
computer defines one or more neighborhoods 244, 246 that include a predetermined number of 
voxels surrounding the pair (step 252). The computer then projects the neighborhoods onto the 
next 2D slice at level N+l by identifying the voxels on the next slice that are immediately 
adjacent the voxels in the neighborhoods on the original slice (step 254). The neighborhoods are 
made large enough to ensure that they include the cusp voxels on the N+l slice. In the example 
of FIG. 13, each cusp voxel 240, 242 lies at the center of a neighborhood 244, 246 of twenty-five 
voxels arranged in a 5 x 5 square. 

In searching for the cusp voxels on the N+l slice, the computer tests the image values for 
all voxels in the projected neighborhoods to identify those associated with the background image 
and those associated with the dentition (step 256). In the illustrated example, voxels in the 
background are black and voxels in the dentition are white. The computer identifies the cusp 
voxels on the N+l slice by locating the pair of black voxels in the two neighborhoods that lie 
closest together (step 258). The computer then repeats this process for all remaining slices (step 
259). 

FIGS. 15 and 16 illustrate another technique, known as "arch curve fitting," for 
identifying interproximal margins between teeth in the dentition. The arch curve fitting 
technique, which also applies to 2D cross-sectional slices of the dentition, involves the creation 
of a curve 260 that fits among the voxels on the 2D cross-sectional surface 262 of the dentition 
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arch 264. A series of closely spaced line segments 268, each bounded by the cross-sectional 
surface 268, are formed along the curve 260, roughly perpendicular to the curve 260, throughout 
the 2D cross section 264. In general, the shortest of these line segments 268 lie on or near the 
interproximal margins; thus computer identifies the cusps that define the interproximal margins 

5 by determining the relative lengths of the line segments 268. 

When applying the arch curve fitting technique, the computer begins by selecting a 2D 
slice (step 270) and identifying the voxels associated with the surface 262 of the cross-sectional 
arch 264 (step 272). The computer then defines a curve 260 that fits among the voxels on the 
surface 262 of the arch (step 274). The computer creates the curve using any of a variety of 

1 0 techniques, a few of which are discussed below. The computer then creates a series of line 
segments that are roughly perpendicular to the curve and are bounded by the cross-sectional 
surface 262 (step 276). The line segments are approximately evenly spaced with a spacing 
distance that depends upon the required resolution and the acceptable computing time. Greater 
resolution leads to more line segments and thus greater computing time. In general, a spacing on 

1 5 the order of 0.4 mm is sufficient in the initial pass of the arch curve fitting technique. 

The computer calculates the length of each line segment (step 278) and then identifies 
those line segments that form local minima in length (step 280). These line segments roughly 
approximate the locations of the interproximal boundaries, and the computer labels the voxels 
that bound these segments as cusp voxels (step 282). The computer repeats this process for each 

20 of the 2D slices (step 284) and then uses the cusp voxels to define 3D cutting surfaces that 
approximate the interproximal margins. 

In some implementations, the computer refines the arch cusp determination by creating 
several additional sets of line segments, each centered around the arch cusps identified on the 
first pass. The line segments are spaced more narrowly on this pass to provide greater resolution 

25 in identifying the actual positions of the arch cusps. 

The computer uses any of a variety of curve fitting techniques to create the curve through 
the arch. One technique involves the creation of a catenary curve with endpoints lying at the two 
ends 265, 267 (FIG. 15) of the arch. The catenary curve is defined by the equation y = a + 
bcosh(cx), and the computer fits the curve to the arch by selecting appropriate values for the 

30 constants a, b, and c. Another technique involves the creation of two curves, one fitted among 
voxels lying on the front surface 271 of the arch, and the other fitted among voxels on the rear 
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surface 273. A third curve, which guides the placement of the line segments above, passes 
through the middle of the area lying between the first two curves. 

FIGS. 17 and 18 illustrate another technique for constructing a curve through the arch. 
This technique involves the creation of a series of initial line segments through the arch 264 and 
5 the subsequent formation of a curve 290 fitted among the midpoints of these line segments This 
curve 290 serves as the arch curve in the arch curve fitting technique described above. 

In applying this technique, the computer first locates an end 265 of the arch (step 300) 
and creates a line segment 291 that passes through the arch 264 near this end 265 (step 301). 
The line segment 291 is bounded by voxels 292a-b lying on the surface of the arch. The 

10 computer then determines the midpoint 293 of the line segment 291 (step 302), selects a voxel 
294 located particular distance from the midpoint 293 (step 304), and creates a second line 
segment 295 that is parallel to the initial line segment 291 and that includes the selected voxel 
294 (step 306). The computer then calculates the midpoint 296 of the second segment 295 (step 
308) and rotates the second segment 295 to the orientation 295 1 that gives the segment its 

1 5 minimum possible length (step 309). In some cases, the computer limits the second segment 295 
to a predetermined amount of rotation (e.g., ± 10°). 

The computer then selects a voxel 297 located a particular distance from the midpoint 
296 of the second segment 295 (step 310) and creates a third line segment 298 that is parallel to 
the second line segment 295 and that includes the selected voxel 297 (step 312). The computer 

20 calculates the midpoint 299 of the third segment 298 (step 3 14) and rotates the segment 298 to 
the orientation 298' that gives the segment its shortest possible length (step 316). The computer 
continues adding line segments in this manner until the other end of the cross-sectional arch is 
reached (step 318). The computer then creates a curve that fits among the midpoints of the line 
segments (step 320) and uses this curve in applying the arch fitting technique described above. 

25 FIGS. 19A, 19B and 20 illustrate an alternative technique for creating 3D surfaces that 

approximate the geometries and locations of the interproximal margins in the patient's dentition. 
This technique involves the creation of 2D planes that intersect the 3D dentition model at 
locations that approximate the interproximal margins. In general, the computer defines a series 
of planes, beginning with an initial plane 330 at one end 331 of the arch 332, that are roughly 

30 perpendicular to the occlusal plane of the dentition model ("Vertical" planes). Each plane 
intersects the dentition model to form a 2D cross section 334. If the planes are spaced 
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sufficiently close to each other, the planes with the smallest cross-sectional areas approximate 
the locations of the interproximal margins in the dentition. The computer locates the 
interproximal regions more precisely by rotating each plane about two orthogonal axes 336, 338 
until the plane reaches the orientation that yields the smallest possible cross-sectional area. 

5 In one implementation of this technique, the computer first identifies one end of the arch 

in the dentition model (step 340). The computer then creates a vertical plane 330 through the 
arch near this end (step 342) and identifies the center point 331 of the plane 330 (step 344). The 
computer then selects a voxel located a predetermined distance from the center point (step 345) 
and creates a second plane 333 that is parallel to the initial plane and that includes the selected 

10 voxel (step 346). The computer calculates the midpoint of the second plane (step 348) and 

rotates the second plane about two orthogonal axes that intersect at the midpoint (step 350). The 
computer stops rotating the plane upon finding the orientation that yields the minimum cross- 
sectional area. In some cases, the computer limits the plane to a predetermined amount of 
rotation (e.g., ± 10° about each axis). The computer then selects a voxel located a particular 

15 distance from the midpoint of the second plane (step 352) and creates a third plane that is parallel 
to the second plane and that includes the selected voxel (step 354). The computer calculates the 
midpoint of the third plane (step 356) and rotates the plane to the orientation that yields the 
smallest possible cross-sectional area (step 357). The computer continues adding and rotating 
planes in this manner until the other end of the arch is reached (step 358). The computer 

20 identifies the planes at which local minima in cross-sectional area occur and labels these planes 
as "interproximal planes," which approximate the locations of the interproximal margins (step 
360). 

One variation of this technique, described in FIG. 2 1 , allows the computer to refine its 
identification of interproximal planes by creating additional, more closely positioned planes in 

25 areas around the planes labeled as interproximal. The computer first creates a curve that fits 
among the midpoints of the planes labeled as interproximal planes (step 372) and then creates a 
set of additional planes along this curve (step 374). The additional planes are not evenly spaced 
along the curve, but rather are concentrated around the interproximal margins. The planes in 
each interproximal area are spaced very closely (e.g., 0.05 mm from each other). The computer 

30 rotates each of the newly constructed planes about two orthogonal axes until the plane reaches its 
minimum cross-sectional area (step 376). The computer then selects the plane in each cluster 
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with the smallest cross-sectional area as the plane that most closely approximates the 
interproximal margin (step 378). 

FIGS. 22, 23, and 24 illustrate a technique for identifying the gingival margin that defines 
the boundary between tooth and gum in the patient's dentition. This technique involves the 
5 creation of a series of vertical 2D planes 380, or slices, that intersect the dentition model roughly 
perpendicular to the occlusal plane (see FIG. 19A). The cross-sectional surface 382 of the 
dentition model in each of these planes 380 includes cusps 384, 386 that represent the gingival 
margin. The computer identifies the gingival margin by applying one or more of the cusp 
detection techniques described above. 

10 One technique is very similar to the neighborhood filtered cusp detection technique 

described above, in that voxel neighborhoods 388, 390 are defined on one of the 2D planes to 
focus the computer's search for cusps on an adjacent 2D plane. Upon detecting a pair of cusps 
384, 386 on one 2D plane (step 400), the computer defines one or more neighborhoods 388, 390 
to include a predetermined number of voxels surrounding the pair (step 402). The computer 

15 projects the neighborhoods onto an adjacent 2D plane by identifying the voxels on the adjacent 
plane that correspond to the voxels in the neighborhoods 388, 390 on the original plane (step 
404). The computer then identifies the pair of black voxels that lie closest together in the two 
neighborhoods on the adjacent plane, labeling these voxels as lying in the cusp (step 406). The 
computer repeats this process for all remaining planes (step 408). 

20 Many of these automated segmentation techniques are even more useful and efficient 

when used in conjunction with human-assisted techniques. For example, techniques that rely on 
the identification of the interproximal or gingival margins function more quickly and effectively 
when a human user first highlights the interproximal or gingival cusps in an image of the 
dentition model. One technique for receiving this type of information from the user is by 

25 displaying a 2D or 3D representation and allowing the user to highlight individual voxels in the 
display. Another technique allows the user to scroll through a series of 2D cross-sectional slices, 
identifying those voxels that represent key features such as interproximal or gingival cusps, as in 
the neighborhood-filtered seed cusp detection technique described above (FIGS. 32, 33, and 34). 
Some of these techniques rely on user interface tools such as cursors and bounding-box markers. 

30 FIGS. 35A-35F illustrate another technique for separating teeth from gingival tissue in 

the dentition model. This technique is a human-assisted technique in which the computer 
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displays an image of the dentition model (step 760) and allows a human operator to identify, for 
each tooth, the gingival margin, or gum line 600, encircling the tooth crown 602 (step 762). 
Some applications of this technique involve displaying a 3D volumetric image of the dentition 
model and allowing the user to select, with an input device such as a mouse, the voxels that 

5 define the gingival line 600 around each tooth crown 602. The computer then uses the identified 
gingival line to model the tooth roots and to create a cutting surface that separates the tooth, 
including the root model, from the gingival tissue 604. 

Once the human operator has identified the gingival line 600, the computer selects a point 
: 606 that lies at or near the center of the tooth crown 602 (step 764). One way of choosing this 

1 0 point is by selecting a 2D image slice that is parallel to the dentition's occlusal plane and that 
intersects the tooth crown 602, and then averaging the x~ and y-coordinate values of all voxels in 
this 2D slice that lie on the surface 608 of the tooth crown 602. After selecting the center point 
606, the computer defines several points 605 on the gingival line 600 (step 766) and fits a plane 
610 among these points 605 (step 768). The computer then creates a line segment 612 that is 

1 5 normal to the plane 610 and that extends a predetermined distance from the selected center point 
606 (step 770). The expected size of a typical tooth or the actual size of the patient's tooth 
determines the length of the line segment 612. A length on the order of two centimeters is 
sufficient to model most tooth roots. The computer defines a sphere 614, or a partial sphere, 
centered at the selected center point 606 (step 772). The radius of the sphere 614 is determined 

20 by the length of the line segment 612. 

The computer then shifts the plane 610 along the line segment 612 so that the plane 610 
is tangential to the sphere 614 (step 774). In some applications, the computer allows the human 
operator to slide the plane 610 along the surface of the sphere 614 to adjust the orientation of the 
plane 610 (step 776). This is useful, for example, when the tooth crown 602 is tilted, which 

25 suggests that the tooth roots also are tilted. The computer then creates a projection 616 of the 
gingival line 600 on the shifted plane 610 (step 778). The tooth roots are modeled by creating a 
surface 618 that connects the gingival line 600 to the projection 616 (step 780). The computer 
uses this surface as a cutting surface to separate the tooth from the gingival tissue. The cutting 
surface extends in a direction that is roughly perpendicular to the occlusal surface of the tooth 

30 crown 602. 
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In general, the surface 618 that connects the gingival line 600 to the projection is formed 
by straight line segments that extend between the gingival line and the projection. However, 
some implementations allow curvature along these line segments. In some applications, the 
computer scales the projection 616 to be larger or smaller than the gingival line 600, which gives 
5 the surface 618 a tapered shape (step 782). Many of these applications allow the computer, with 
or without human assistance, to change the profile of the tapered surface so that the rate of 
tapering changes along the length of the surface 618 (step 784). For example, some surfaces 
taper more rapidly as distance from the tooth crown increases. 

FIGS. 37A-C and 38 illustrate another human-assisted technique for separating teeth 

10 from gingival tissue in the dentition model. This technique involves displaying an image of the 
dentition model to a human operator (step 790) and allowing the operator to trace the gingival 
lines 620, 622 on the buccal and lingual sides of the dental arch (step 792). This produces two 
3D curves 624, 626 representing the gingival lines 620, 622 on the buccal and lingual surfaces. 
The computer uses these curves 624, 626 to create a 3D cutting surface 628 that separates the 

15 tooth crowns 630, 632 from the gingival tissue 634 in the dentition model (step 794). The 
cutting surface 628 is roughly parallel to the occlusal surface of the tooth crowns 630, 632. 

FIGS. 37C and 39 illustrate one technique for defining the cutting surface 628. In 
general, the computer creates the cutting surface 628 by defining points 636, 638 along each of 
the 3D curves 624, 626 and defining the cutting surface 628 to fit among the points 636, 638. 

20 The computer first defines the points 636, 638 on the 3D curves 624, 626 (step 800) and then 
defines a point 640 at or near the center of each tooth crown 630 (step 802). One way of 
defining the center point 640 is by averaging the x-, y-, and z-coordinate values for all of the 
points 636, 638 lying on the portions of the gingival curves 624, 626 associated with that tooth. 
The computer then creates a triangular surface mesh 642 using the center point 640 and the 

25 points 636, 638 on the gingival curves as vertices (step 804). The computer uses this surface 
mesh 642 to cut the tooth crowns away from the gingival tissue (step 806). In some 
implementations, a tooth root model is created for each crown, eg., by projecting the gingival 
curves onto a distant plane, as described above (step 808). The computer connects the roots to 
the crowns to complete the individual tooth models (step 810). 

30 All of the segmentation techniques described above are useful in creating digital models 

of individual teeth, as well as a model of gingival tissue surrounding the teeth. In some cases, 
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the computer identifies and segments the teeth using one of these techniques to form the 
individual tooth models, and then uses all remaining data to create the gingival model. 

Other Implementations 

5 In many instances, the computer creates proposals for segmenting the dentition model 

and then allows the user to select the best alternative. For example, one version of the arch curve 
fitting technique described above requires the computer to create a candidate catenary or spline 
curve, which the user is allowed to modify by manipulating the mathematical control parameters. 
Other techniques involve displaying several surfaces that are candidate cutting surfaces and 

10 allowing the user to select the appropriate surfaces. 

Some implementations of the invention are realized in digital electronic circuitry, such as 
an application specific integrated circuit (ASIC); others are realized in computer hardware, 
firmware, and software, or in combinations of digital circuitry and computer components. The 
invention is usually embodied, at least in part, as a computer program tangibly stored in a 

1 5 machine-readable storage device for execution by a computer processor. In these situations, 
methods embodying the invention are performed when the processor executes instructions 
organized into program modules, operating on input data and generating output Suitable 
processors include general and special purpose microprocessors, which generally receive 
instructions and data from read-only memory and/or random access memory devices. Storage 

20 devices that are suitable for tangibly embodying computer program instructions include all forms 
of nonvolatile memory, including semiconductor memory devices, such as EPROM, EEPROM, 
and flash memory devices; magnetic disks such as internal hard disks and removable disks; 
magneto-optical disks; and CD-ROM. 

The invention has been described in terms of particular embodiments. Other 

25 embodiments are within the scope of the following claims. 
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CLAIMS 

1 . A computer-implemented method for use in creating a digital model of an 
individual component of a patient's dentition, the method comprising: 

5 (a) receiving a data set that forms a three-dimensional (3D) representation of 

the patient's dentition; 

(b) applying a computer-implemented test to the data set to identify data 
elements that represent portions of an individual component of the patient's dentition; 
and 

1 0 (c) creating a digital model of the individual component based upon the 

identified data elements. 

2. The method of claim 1 , wherein the data set includes data taken from at least one 
of the following sources: two-dimensional (2D) x-ray data and three-dimensional (3D) x-ray 

15 data. 

3. The method of claim 1 , wherein the data set includes data taken from at least one 
of the following sources: computed tomography (CT) scan data and magnetic resonance imaging 
(MRI) scan data. 

20 

4. The method of claim 1, wherein the data set includes data taken from a 
photographic image of the patient's dentition. 

5. The method of claim 1 , wherein some of the data is obtained by imaging a 
25 physical model of the patient' s teeth. 

6. The method of claim 1, wherein some of the data is obtained by imaging the 
patient's teeth directly. 

30 7. The method of claim 1 , wherein the data set forms a 3D volumetric representation 

of the patient's dentition. 
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8. The method of claim 1 , wherein the data set includes geometric surface data that 
forms a 3D geometric surface model of the patient's dentition. 

9. The method of claim 1, wherein the individual component is an individual tooth 
in the patient' s dentition. 

10. The method of claim 1, wherein the individual component includes gum tissue 
found in the patient's dentition. 

1 1 . The method of claim 1 , wherein applying the computer-implemented test includes 
receiving information input by a human user to identify a boundary of the individual component 
to be modeled. 

12. The method of claim 11, wherein receiving information includes receiving 
position data from a computer-implemented tool through which the human user identifies the 
boundary in a graphical representation of the patient's dentition. 

13. The method of claim 12, wherein the computer-implemented tool is a saw tool 
that allows the user to identify the boundary by defining a curve in the graphical representation 
that separates the data elements associated with the individual component from other elements of 
the data set 

14. The method of claim 12, wherein the computer-implemented tool is an eraser tool 
that allows the user to identify the boundary by erasing a portion of the graphical representation 
representing the boundary. 

15. The method of claim 1 , wherein receiving the data, applying the computer- 
implemented test, and creating the electronic model all are carried out by a computer without 
human intervention. 
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1 6. The method of claim 1 , wherein applying the computer-implemented test includes 
automatically applying a rule to identify a boundary of the individual component to be modeled. 

1 7. The method of claim 1 6, wherein the boundary includes a surface of a tooth. 

5 

1 8. The method of claim 1 6, wherein the boundary includes a gingival margin. 

19. The method of claim 1 , wherein applying the computer-implemented test includes 
identifying elements of the data set that represent a structural core of the individual component to 

1 0 be modeled and labeling those data elements as belonging to the individual component 

20. The method of claim 1 9, wherein the individual component to be modeled 
includes an individual tooth and the structural core approximately coincides with neurological 
roots of the tooth. 

15 

2 1 . The method of claim 1 9, wherein applying the computer-implemented test 
includes applying a test to link other data elements to those representing the structural core and 
labeling the linked data elements as belonging to the individual component 

20 22. The method of claim 21, wherein applying the test to link other data elements to 

those representing the structural core includes assigning a distance measure to each element of 
the data set, where the distance measure indicates a measured distance between a reference point 
in the dentition and the portion of the dentition represented by the data element to which the 
distance measure is assigned. 

25 

23. The method of claim 22, wherein applying the test to link other data elements 
includes linking a data element to the structural core if the assigned distance measure is less than 
the distance measure assigned to a data element representing a portion of the structural core. 

30 24. The method of claim 22, wherein the reference point lies on a tooth surface. 
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25. The method of claim 21, wherein applying the test to link other data elements to 
the structural core includes applying a test to determine whether a data element lies outside of the 
dentition and, if so, labeling the data element as a background element. 

5 26. The method of claim 25, wherein applying the test to determine whether the data 

element lies outside of the dentition includes comparing an image value associated with the data, 
element to a threshold value. 



27. The method of claim 19, further comprising applying another computer- 

10 implemented test to identify dements of the data set that represent a structural core of another 
individual component of the dentition and labeling those data elements as belonging to the other 
individual component. 

28. The method of claim 27, wherein applying the computer-implemented tests 

1 5 includes applying tests to link other elements of the data set to those representing the structural 
cores of the individual components and labeling the linked elements as belonging to the 
individual components to which they are linked. 

29. The method of claim 28, wherein applying the tests to link other data elements to 
20 the structural cores of the individual components includes determining whether a data element 

already is labeled as belonging to one of the individual components. 

30. The method of claim 1, wherein applying the computer implemented test includes 
identifying an initial 2D cross-section of the individual component having continuous latitudinal 

25 width, a relative minimum value of which occurs at an end of the initial cross-section. 

31. . . The method of claim 30, wherein applying the computer-implemented test 
includes isolating portions of the data corresponding to the initial 2D cross-section of the 
individual component to be modeled. 

30 
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32. The method of claim 3 1 , wherein the received data includes 3D image data 
obtained by imaging the individual component volumetrically, and wherein isolating portions of 
the data corresponding to the initial 2D cross-section includes isolating elements of the 3D image 
data representing the initial 2D cross-section. 

5 ' 

33. The method of claim 30, wherein applying the computer-implemented test 
includes applying a test to identify the end of the initial cross-section at which the relative 
minimum value of the latitudinal width occurs. 

10 34. The method of claim 33, wherein applying the test to identify the end of the initial 

cross-section includes: 

(a) establishing line segments within the initial cross-section, each of which is 
bounded at each end by an endpoint lying on a surface of the individual component, and 
each of which is roughly perpendicular to a latitudinal axis of the individual component; 
15 (b) calculating a length for each line segment; and 

(c) identifying elements of the data set that correspond to the endpoints of the 
line segment with the shortest length. 

35. The method of claim 34, wherein applying the computer-implemented test also 
20 includes: 

(a) isolating portions of the data set corresponding to other 2D cross-sections 
of the individual component, all lying in planes parallel to the initial 2D cross-section; 

(b) for each of the other cross-sections, identifying data elements that 
correspond to endpoints of a line segment representing an end of the cross-section; and 

25 (c) defining a solid surface that contains all of the identified data elements. 

36. The method of claim 35, further comprising labeling the solid surface as 
representing a surface of the individual component to be modeled. 

30 37. The method of claim 35, further comprising using the data elements identified in 

the initial cross-section as guides for identifying the data elements in the other cross-sections. 



28 



WO 00/19935 PCT/US99/23532 



' 38. The method of claim 34, wherein applying the test to identify the end of the initial 
cross-section includes first creating an initial curve that is roughly perpendicular to the latitudinal 
axis of the individual component and that is fitted between the surfaces of the 2D cross-section 
5 on which the endpoints of the line segments will lie. 

39. The method of claim 38, wherein establishing the line segments includes first 
establishing a set of initial line segments that are roughly perpendicular to the curve and to the 
latitudinal axis and that have endpoints lying on the surfaces of the individual component 

10 

40. The method of claim 39, wherein establishing the line segments also includes 
pivoting each initial line segment about a point at which the initial line segment intersects the 
curve until the initial line segment has its shortest possible length. 

15 41 . The method of claim 40, wherein establishing the line segments also includes: 

(a) locating a midpoint for each of the initial line segments after pivoting; and 

(b) creating a refined curve that passes through all of the midpoints. 

42. The method of claim 41, wherein establishing the line segments also includes 
20 creating the line segments to be perpendicular to the refined curve. 

43. The method of claim 38j wherein the individual component is a tooth and the 
curve is a portion of a larger curve fitted among the lingual and buccal surfaces of all teeth in a 
2D cross-section of a tooth arch in which the tooth lies. 



25 



44. The method of claim 43, wherein the larger curve is a catenary. 

45. The method of claim 43, wherein the larger curve is created by manipulating 
mathematical control points to fit the curve to the shape of the cross-section of the tooth arch. 



30 
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46. The method of claim 34, wherein establishing the line segments includes first 
establishing an initial line segment by creating a line that intersects the initial 2D cross-section, 
such that the initial line segment has endpoints that lie on surfaces of the individual component 

47. The method of claim 46, wherein establishing the line segments also includes 
establishing at least one additional line segment parallel to and spaced a predetermined distance 
from a previously established line segment 

48. The method of claim 47, wherein establishing the line segments also includes, for 
each additional line segment, locating a midpoint of the additional line segment and pivoting the 
additional line segment about the midpoint until the additional line segment has its shortest 
possible length. 

49. The method of claim 48, wherein establishing the line segments also includes 
limiting the rotation of each additional line segment to no more than a predetermined amount 

50. The method of claim 49, wherein the rotation of each additional line segment is 
limited to no more than approximately +/- 10°. 

51 . The method of claim 48, wherein establishing the line segments also includes 
establishing a curve that is fitted among the midpoints of the additional line segments. 

52. The method of claim 51, wherein establishing the line segments includes 
establishing the line segments to be perpendicular to the curve. 

53. The method of claim 52, wherein establishing the line segments includes locating 
midpoints for each of the line segments and pivoting each line segment about its midpoint until 
the line segment has its shortest possible length. 
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54. The method of claim 30, wherein the individual component is a tooth and the 
relative minimum value of the initial 2D cross-section lies on an interproximal surface of the 
tooth. 

5 55. The method of claim 54, wherein identifying the initial 2D cross-section includes 

isolating elements of the data set that correspond to 2D cross-sections of the tooth lying in 
parallel planes between the roots and the occlusal surface of the tooth. 

56. The method of claim 55, wherein identifying the initial 2D cross-section also 

10 includes identifying adjacent ones of the 2D cross-sections in which the interproximal surface of 
the tooth is obscured by gum tissue in one of the adjacent cross-sections and is not obscured by 
gum tissue in the other adjacent cross-section. 

57. The method of claim 56, wherein identifying the initial 2D cross-section also 
15 includes selecting as the initial 2D cross-section the adjacent cross-section in which the 

interproximal surface of the tooth is not obscured by gum tissue. 

58. The method of claim 55, wherein identifying the initial 2D cross-section also 
includes, for each of the isolated cross-sections, establishing a contour line that outlines the 

20 shape of the dentition in that cross-section. 

59. The method of claim 58, wherein identifying the initial 2D cross-section also, 
includes applying a test to each of the isolated cross-sections to identify those cross-sections in 
which the interproximal surface of the tooth is not obscured by gum tissue. 

25 

60. The method of claim 59, wherein applying the test includes calculating the rate of 
curvature of the contour line. 

6 1 . The method of claim 59, wherein identifying the initial 2D cross-section includes 
30 selecting as the initial 2D cross-section the isolated cross-section that lies closest to the roots of 

the tooth and in which the interproximal surface of the tooth is not obscured by gum tissue. 
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62. The method of claim 30, wherein applying the computer-implemented test also 
includes identifying two elements of the data set that define endpoints of a line segment spanning 
the relative minimum width of the initial 2D cross-section. 

5 

63. The method of claim 62, wherein applying the computer-implemented test also 
includes defining, for each endpoint, a neighborhood containing a predetermined number of 
elements of the data set near the endpoint in the initial 2D cross-section. 

1 0 64. The method of claim 63, wherein applying the computer-implemented test also 

includes identifying an additional 2D cross-section of the individual component in a plane 
parallel and adjacent to the initial 2D cross-section, where the additional 2D cross-section also 
has a continuous, latitudinal width with a relative minimum value occurring at one end of the 
cross-section. 

15 

65. The method of claim 64, wherein applying the computer-implemented test also 
includes identifying two elements of the data set that define endpoints of a line segment spanning 
the relative minimum width of the additional 2D cross-section by: 

(a) defining two neighborhoods of data elements, each containing elements of 
20 the data set that are adjacent to the data elements contained in the neighborhoods defined 

for the initial 2D cross-section; and 

(b) identifying one data element in each neighborhood that corresponds to one 
of the endpoints of the line segment spanning the relative minimum width of the 
additional 2D cross-section. 

25 

66. The method of claim 65, further comprising establishing a solid surface that is 
fitted among line segments spanning the relative minimum widths of the parallel 2D cross- 
sections. 

30 67. The method of claim 66, wherein the individual component to be modeled is a 

tooth and the solid surface represents an interproximal surface of the tooth. 
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68. The method of claim 30, further comprising receiving information provided by a 
human user that identifies elements of the data set that correspond to the relative minimum width 
of the initial 2D cross-section. 

5 

69. The method of claim 68, further comprising displaying a graphical representation 
of the patient's dentition in which the user identifies portions corresponding to the relative 
minimum width of the cross-section. 

10 70. The method of claim 69, wherein the graphical representation is three 

dimensional. 

71 . The method of claim 69, wherein the graphical representation includes a 2D 
representation of the initial 2D cross-section. 

15 

72. The method of claim 71, further comprising receiving the information from an 
input device used by the human user to identify the relative minimum width of the initial 2D 
cross-section in the graphical representation. 

20 73. The method of claim 71, wherein the initial 2D cross-section is one of many 2D 

cross-sections displayed to the human user. 

74. The method of claim 7 1 , further comprising receiving information from the 
human user identifying which of the displayed 2D cross-sections is the initial 2D cross-section. 

25 

75. A computer-implemented method for use in creating a digital model of a tooth in 
a patient's dentition, the method comprising: 

(a) receiving a three-dimensional (3D) data set representing the patient's 
dentition; 

30 (b) applying a computer-implemented test to identify data elements that 

represent an interproximal margin between two teeth in the dentition; 
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(c) applying another computer-implemented test to select data elements that 
lie on one side of the interproximal margin for inclusion in the digital model. 

76. The method of claim 75, further comprising creating a set of 2D planes that 

5 intersect the dentition roughly perpendicular to an occlusal plane of the dentition, each 2D plane 
including data elements that form a 2D cross-section of the dentition. 

77. The method of claim 76, further comprising identifying the 2D plane with the 
smallest cross-sectional area. 

10 

78. The method of claim 77, further comprising rotating the 2D plane with the 
smallest cross-sectional area to at least one other orientation to form at least one other 2D cross- 
section of the dentition. 

15 79. The method of claim 78, further comprising selecting the orientation that gives 

the rotated plane its smallest possible cross-sectional area. 

80. The method of claim 79, further comprising identifying data elements that 
represent the selected orientation of the rotated plane as lying on an interproximal margin. 

20 

81 . The method of claim 78, wherein the plane is rotated about two orthogonal lines 
passing through its center point. 

82. The method of claim 77, further comprising creating a set of additional 2D planes 
25 in the vicinity of the 2D plane with the smallest cross-sectional area. 

83. The method of claim 82, further comprising identifying the plane in the set of 
additional planes that has the smallest cross-sectional area. 
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84. The method of claim 83, further comprising rotating the plane with the smallest 
cross-sectional area to at least one other orientation to form at least one other 2D cross-section of 
the dentition. 

85. The method of claim 84, further comprising selecting the orientation that 
produces the 2D cross-section with the smallest possible area. 

86. The method of claim 76, wherein creating a set of 2D planes includes creating an 
initial plane near one end of the dentition. 

87. The method of claim 86, further comprising selecting a point in the dentition that 
is a predetermined distance from the initial plane and creating a second plane. 

88. The method of claim 87, wherein the second plane is roughly parallel to the initial 

plane. 

89. The method of claim 87, further comprising rotating the second plane to at least 
one additional orientation to form at least one additional 2D cross-section of the dentition. 

90. The method of claim 89, further comprising selecting the orientation that 
produces the 2D cross-section with the smallest cross-sectional area. 

91 . The method of claim 89, further comprising selecting a point that is a 
predetermined distance from the second plane and creating a third plane that includes the 
selected point. 

92. The method of claim 91, further comprising rotating the third plane to at least one 
other orientation to create at least one additional 2D cross-section of the dentition. 
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93. The method of claim 91, further comprising creating additional planes, each 
including a point that is a predetermined distance from a preceding plane, until the other end of 
the dentition is reached. 

94. The method of claim 93, further comprising identifying at least one plane having 
a local minimum in cross-sectional area. 

95. The method of claim 93, further comprising identifying a centerpoint of the cross- 
section in each of the planes and creating a curve that fits among the identified centeipoints. 

96. The method of claim 95, further comprising creating a set of additional 2D planes 
along the curve, where the curve is roughly normal to each of the additional planes, and where 
each of the additional planes is roughly perpendicular to the occlusal plane. 

97. The method of claim 96, further comprising identifying at least one of the 
additional planes that has a local minimum in cross-sectional area. 

98. A computer-implemented method for use in creating a digital model of a tooth in 
a patient's dentition, the method comprising: 

(a) receiving a 3D dataset representing at least a portion of the patient's 
dentition, including at least a portion of a tooth and gum tissue surrounding the tooth; 

(b) applying a test to identify data elements lying on a gingival boundary that 
occurs where the tooth and the gum tissue meet; and 

(c) applying a test to the data elements lying on the boundary to identify other 
data elements representing portions of the tooth. 

99. The method of claim 98, wherein applying the test to identify data elements on 
the gingival boundary includes creating an initial 2D plane that intersects the dentition roughly 
perpendicular to an occlusal plane of the dentition and that includes data elements representing 
an initial cross-sectional surface of the dentition. 
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100. The method of claim 99, wherein applying the test includes locating a cusp in the 
initial cross-sectional surface. 

101. The method of claim 100, wherein locating the cusp includes calculating rate of 
5 curvature of the initial cross-sectional area at selected points on the cross-sectional surface. 

102. The method of claim 101, wherein locating the cusp includes identifying the point 
at which the rate of curvature is greatest. 

10 : 103. The method of claim 100, wherein applying the test includes creating a second 2D 

plane that is roughly parallel to the initial 2D plane and that includes data elements representing 
a second cross-sectional surface of the dentition. 

1 04. The method of claim 1 03, wherein applying the test includes locating a cusp in 
15 the second cross-sectional surface. 

105. The method of claim 104, wherein locating the cusp in the second cross-sectional 
surface includes defining a neighborhood of data elements around the cusp in the initial cross- 
sectional surface and projecting the neighborhood onto the second cross-sectional surface. 

20 

106. The method of claim 105, wherein locating the cusp in the second cross-sectional 
surface includes searching for the cusp only within the neighborhood projected onto the second 
cross-sectional surface. 

25 1 07. The method of claim 99, wherein applying the test includes locating two cusps in 

the initial cross-sectional surface. 



1 08. The method of claim 107, wherein applying the test includes creating a second 2D 
plane that is roughly parallel to the initial 2D plane and that includes data elements representing 
30 a second cross-sectional surface of the dentition. 
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109. The method of claim 108, wherein applying the test includes locating two cusps in 
the second cross-sectional surface. 

110. The method of claim 1 09, wherein locating the cusps in the second cross-sectional 
5 surface includes defining two neighborhoods of data elements around the two cusps in the initial 

cross-sectional surface and projecting the neighborhoods onto the second cross-sectional surface. 

111. The method of claim 110, wherein each neighborhood projected onto the second 
cross-sectional surface includes data elements representing portions of the tooth and data 

10 elements representing the gum tissue surrounding the tooth. 

1 12. The method of claim 111, wherein the data elements representing the tooth 
include voxels of one color and the data elements representing the gum tissue include voxels of 
another color. 

15 

1 13. The method of claim 111, wherein locating the cusps in the second cross-sectional 
surface includes locating the pair of data elements representing gum tissue that lie closest 
together, where each of the two neighborhoods projected onto the second cross-sectional surface 
includes one of the data elements in the pair. 

20 

1 14. The method of claim 98, wherein applying the test to identify data elements on 
the gingival boundary includes creating a series of roughly parallel 2D planes, each intersecting 
the dentition roughly perpendicular to an occlusal plane of the dentition, and each including data 
elements that represent a cross-sectional surface of the dentition. 

25 

115. The method of claim 1 14, wherein the cross-sectional surface in each 2D plane 
includes two cusps that roughly identify the locations of the gingival boundary. 

116. The method of claim 115, wherein applying the test includes identifying the cusps 
30 in each cross-sectional surface. 
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1 17. The method of claim 1 16, wherein identifying the cusps includes locating the 
cusps in one of the planes and then confining the search for cusps in an adjacent plane to a 
predetermined area in the vicinity of the identified cusps. 

5 118. The method of claim 1 1 5, further comprising allowing a human user to select data 

elements that roughly identify the locations of the cusps in a selected one of the cross-sectional . 
areas. 

1 19. The method of claim 1 1 8, further comprising searching for cusps in the selected 
1 0 cross-sectional area and confining the search to data elements lying within a predetermined area 

in the vicinity of the data elements selected by the human user. 

120. The method of claim 1 1 8, further comprising searching for cusps in an adjacent 
cross-sectional area and confining the search to data elements lying within a predetermined area 

15 in the vicinity of the data elements selected by the human user. 

121. A computer program, stored on a tangible storage medium, for use in creating a 
digital model of an individual component of a patient's dentition, the program including 
executable instructions that, when executed by a computer, cause the computer to: 

20 (a) receive a data set that forms a three-dimensional (3D) representation of the 

patient's dentition; 

(b) apply a test to the data set to identify data elements that represent portions 
of an individual component of the patient's dentition; and 

(c) create a digital model of the individual component based upon the 
25 identified data elements. 

122. The program of claim 121 , wherein the computer receives the data, applies the 
test, and creates the electronic model without human intervention. 

30 123. The program of claim 121, wherein the computer, in applying the test, applies a 

rule to identify a boundary of the individual component to be modeled. 
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124. The program of claim 121, wherein the computer, in applying the test, identifies 
elements of the data set that represent a structural core of the individual component to be 
modeled and labels those data elements as belonging to the individual component. 

125. The program of claim 124, wherein the computer, in applying the test, links other 
data elements to those representing the structural core and labels the linked data elements as 
belonging to the individual component. 

126. The program of claim 125, wherein the computer, in applying the test: 

(a) assigns a distance measure to each element of the data set, where the 
distance measure indicates a measured distance between a reference point in the dentition 
and the portion of the dentition represented by the data element to which the distance 
measure is assigned; and 

(b) links a data element to the structural core if the assigned distance measure 
is less than the distance measure assigned to a data element representing a portion of the 
structural core. 

127. The program of claim 121, wherein the computer, in applying the test, identifies 
20 an initial 2D cross-section of the individual component having continuous latitudinal width, a 

relative minimum value of which occurs at an end of the initial cross-section. 



10 



15 



128. The program of claim 127, wherein the computer, in applying the test, identifies 
the end of the initial cross-section at which the relative minimum value of the latitudinal width 
25 occurs by: 

(a) establishing line segments within the initial cross-section, each of which is 
bounded at each end by an endpoint lying on a surface of the individual component, and 
each of which is roughly perpendicular to a latitudinal axis of the individual component; 

(b) calculating a length for each line segment; and 

30 (c) identifying elements of the data set that correspond to the endpoints of the 

line segment with the shortest length. 
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129. The program of claim 128, wherein the computer, in applying the test 

(a) isolates portions of the data set corresponding to other 2D cross-sections 
of the individual component, all lying in planes parallel to the initial 2D cross-section; 

5 (b) for each of the other cross-sections, identifies data elements that 

correspond to endpoints of a line segment representing an end of the cross-section; and . 
(c) defines a solid surface that contains all of the identified data elements. 

130. The program of claim 128, wherein the computer, in applying the test: 

10 (a) first creates an initial curve that is roughly perpendicular to the latitudinal 

axis of the individual component and that is fitted between the surfaces of the 2D cross- 
section on which the endpoints of the line segments lie; 

(b) establishes a set of initial line segments that are roughly perpendicular to 
. the curve and to the latitudinal axis and that have endpoints lying on the surfaces of the 

1 5 individual component; 

(c) pivots each initial line segment about a point at which the initial line 
segment intersects the curve until the initial line segment has its shortest possible length; 

(d) locates a midpoint for each of the initial line segments after pivoting; and 

(e) creates a refined curve that passes through all of the midpoints and that is 
20 roughly normal to all of the line segments. 

131. The program of claim 1 28, wherein the computer, in applying the test, also: 

(a) establishes an initial line segment by creating a line that intersects the 
initial 2D cross-section, such that the initial line segment is bounded by endpoints that lie 

25 on surfaces of the individual component; 

(b) establishes at least one additional line segment parallel to and spaced a 
predetermined distance from a previously established line segment; and 

(c) for each additional line segment, locates a midpoint of the additional line 
segment and pivots the additional line segment about the midpoint until the additional 

30 line segment has its shortest possible length. 



41 



WO 00/19935 PCT/US99/23532 



132. The program of claim 131, wherein the computer, in applying the test, also: 

(a) establishes a curve that is fitted among the midpoints of the additional line 
segments; 

(b) establishes a new set of line segments that are perpendicular to the curve; 

(c) locates midpoints for each of the line segments in the new set; and 

(d) pivots each line segment in the new set about its midpoint until the line 
segment has its shortest possible length. 



133. The program of claim 127, wherein the individual component is a tooth and the 
10 relative minimum value of the initial 2D cross-section lies on an interproximal surface of the 

tooth. 

134. The program of claim 133, wherein the computer, in identifying the initial 2D 
cross-section, isolates elements of the data set that correspond to 2D cross-sections of the tooth 

1 5 lying in parallel planes between the roots and the occlusal surface of the tooth. 



135. The program of claim 134, wherein the computer, in identifying the initial 2D 
cross-section, identifies adjacent ones of the 2D cross-sections in which the interproximal 
surface of the tooth is obscured by gum tissue in one of the adjacent cross-sections and is not 

20 obscured by gum tissue in the other adjacent cross-section. 

136. The program of claim 135, wherein the computer, in identifying the initial 2D 
cross-section, selects as the initial 2D cross-section the adjacent cross-section in which the 
interproximal surface of the tooth is not obscured by gum tissue. 

25 

137. The program of claim 1 34, wherein the computer, in identifying the initial 2D 
cross-section, identifies for each of the isolated cross-sections a contour line that outlines the 
shape of the dentition in that cross-section. 
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138. The program of claim 1 37, wherein the computer, in identifying the initial 2D 
cross-section, applies a test to each of the isolated cross-sections to identify those cross-sections 
in which the interproximal surface of the tooth is not obscured by gum tissue. 

5 139. The program of claim 138, wherein the computer, in applying the test to each of 

the isolated cross-sections, calculates the rate of curvature of the contour line. 

140. The program of claim 138, wherein the computer, in identifying the initial 2D 
cross-section, selects as the initial 2D cross-section the isolated cross-section that lies closest to 

10 the roots of the tooth and in which the interproximal surface of the tooth is not obscured by gum 
tissue. 

141. The program of claim 127, wherein the computer, in applying the test, identifies 
two elements of the data set that define endpoints of a line segment spanning the relative 

15 minimu m width of the initial 2D cross-section. 

142. The program of claim 141, wherein the computer, in applying the test, defines for 
each endpoint a neighborhood containing a predetermined number of elements of the data set 
near the endpoint in the initial 2D cross-section. 

20 

143. The program of claim 142, wherein the computer, in applying the test, identifies 
an additional 2D cross-section of the individual component in a plane parallel and adjacent to the 
initial 2D cross-section, where the additional 2D cross-section also has a continuous, latitudinal 
width with a relative minimum value occurring at one end of the cross-section. 

25 

144. The program of claim 143, wherein the computer, in applying the test, identifies 
two elements of the data set that define endpoints of a line segment spanning the relative 
minimum width of the additional 2D cross-section by: 

(a) defining two neighborhoods of data elements, each containing elements of 
30 the data set that are adjacent to the data elements contained in the neighborhoods defined 

for the initial 2D cross-section; and 
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(b) identifying one data element in each neighborhood that corresponds to one 
of the endpoints of the line segment spanning the relative minimum width of the 
additional 2D cross-section. 

145. The program of claim 144, wherein the computer also establishes a solid surface 
that is fitted among line segments spanning the relative minimum widths of the parallel 2D 
cross-sections. 

146. The program of claim 145, wherein the individual component to be modeled is a 
tooth and the solid surface represents an interproximal surface of the tooth. 

147. A computer program, stored on a tangible storage medium, for use in creating a 
digital model of tooth in a patient's dentition, the program including executable instructions that, 
when executed by a computer, cause the computer to: 

(a) receive a three-dimensional (3D) data set representing the patient's 
dentition; 

(b) apply a test to identify data elements that represent an interproximal 
margin between two teeth in the dentition; 

(c) apply another test to select data elements that lie on one side of the 
interproximal margin for inclusion in the digital model. 

148. The program of claim 147, wherein the computer creates a set of 2D planes that 
intersect the dentition roughly perpendicular to an occlusal plane of the dentition, each 2D plane 
including data elements that form a 2D cross-section of the dentition. 

149. The program of claim 148, wherein the computer identifies the 2D plane with the 
smallest cross-sectional area. 

1 50. The program of claim 1 49, wherein the computer rotates the 2D plane with the 
smallest cross-sectional area to at least one other orientation to form at least one other 2D cross- 
section of the dentition. 
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151. The program of claim 1 50, wherein the computer selects the orientation that gives 
the rotated plane its smallest possible cross-sectional area. 

5 152. The program of claim 151, wherein the computer identifies data elements that , 

represent the selected orientation of the rotated plane as lying on an interproximal margin. 

1 53. The program of claim 1 50, wherein the computer rotates the plane about two 
orthogonal lines passing through its center point. 

10 

1 54. The program of claim 1 49, wherein the computer creates a set of additional 2D 
planes in the vicinity of the 2D plane with the smallest cross-sectional area. 

155. The program of claim 1 54, wherein the computer identifies the plane in the set of 
15 additional planes that has the smallest cross-sectional area. 

1 56. The program of claim 1 55, wherein the computer rotates the plane with the 
smallest cross-sectional area to at least one other orientation to form at least one other 2D cross- 
section of the dentition. 

20 

157. The program of claim 156, wherein the computer selects the orientation that 
produces the 2D cross-section with the smallest possible area. 

158. The program of claim 148, wherein the computer, in creating the set of 2D planes, 
25 creates an initial plane near one end of the dentition. 

1 59. The program of claim 1 58, wherein the computer selects a point in the dentition 
that is a predetermined distance from the initial plane and creates a second plane that includes the 
selected point 

30 
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160. The program of claim 159, wherein the second plane is roughly parallel to the 
initial plane. 

161. The program of claim 1 59, wherein the computer rotates the second plane to at 
5 least one additional orientation to form at least one additional 2D cross-section of the dentition. 

162. The program of claim 161, wherein the computer selects the orientation that 
produces the 2D cross-section with the smallest cross-sectional area. 

10 163. The program of claim 161, wherein the computer selects a point that is a 

predetermined distance from the second plane and creates a third plane that includes the selected 
point. 

164. The program of claim 163, wherein the computer rotates the third plane to at least 
15 one other orientation to create at least one additional 2D cross-section of the dentition. 

1 65 . The program of claim 1 63 , wherein the computer creates additional planes, each 
including a point that is a predetermined distance from a preceding plane, until the other end of 
the dentition is reached. 

20 

166. The program of claim 165, wherein the computer identifies at least one plane 
having a local minimum in cross-sectional area. 

1 67. The program of claim 1 65, wherein the computer identifies a centerpoint of the 
25 cross-section in each of the planes and creates a curve that fits among the identified centerpoints. 

1 68. The program of claim 1 67, wherein the computer creates a set of additional 2D 
planes along the curve, where the curve is roughly normal to each of the additional planes, and 
where each of the additional planes is roughly perpendicular to the occlusal plane. 

30 
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169. The program of claim 168, wherein the computer identifies at least one of the 
additional planes that has a local minimum in cross-sectional area. 

170. A computer program, stored on a tangible storage medium, for use in creating a 
5 digital model of a tooth in a patient's dentition, the program including executable instructions 

that, when executed by a computer, cause the computer to: 

(a) receive a 3D data set representing at least a portion of the patient's 
dentition, including at least a portion of a tooth and gum tissue surrounding the tooth; 

(b) apply a test to identify data elements lying on a gingival boundary that 
10 occurs where the tooth and the gum tissue meet; and 

(c) apply a test to the data elements lying on the boundary to identify other 
data elements representing portions of the tooth. 

1 7 L The program of claim 1 70, wherein the computer, in applying the test to identify 
15 data elements on the gingival boundary, creates an initial 2D plane that intersects the dentition 
roughly perpendicular to an occlusal plane of the dentition and that includes data elements 
representing an initial cross-sectional surface of the dentition. 

172. The program of claim 171, wherein the computer locates a cusp in the initial 
20 cross-sectional surface. 

173. The program of claim 172, wherein the computer, in locating the cusp, calculates 
rate of curvature of the initial cross-sectional area at selected points on the cross-sectional 
surface. 

25 

» 174. The program of claim 173, wherein the computer, in locating the cusp, identifies 
the point at which the rate of curvature is greatest 

175. The program of claim 172, wherein the computer creates a second 2D plane that 
30 is roughly parallel to the initial 2D plane and that includes data elements representing a second 
cross-sectional surface of the dentition. 
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1 76. The program of claim 1 75, wherein the computer locates a cusp in the second 
cross-sectional surface. 

1 77. The program of claim 1 76, wherein the computer, in locating the cusp in the 
second cross-sectional surface, defines a neighborhood of data elements around the cusp in the . 
initial cross-sectional surface and projects the neighborhood onto the second cross-sectional 
surface. 

1 78. The program of claim 1 77, wherein the computer, in locating the cusp in the 
second cross-sectional surface, searches for the cusp only within the neighborhood projected 
onto the second cross-sectional surface. 

1 79. The program of claim 171, wherein the computer locates two cusps in the initial 
cross-sectional surface. 

1 80. The program of claim 1 79, wherein the computer creates a second 2D plane that 
is roughly parallel to the initial 2D plane and that includes data elements representing a second 
cross-sectional surface of the dentition. 

181. The program of claim 180, wherein the computer locates two cusps in the second 
cross-sectional surface. 

1 82. The program of claim 181, wherein the computer, in locating the cusps in the 
second cross-sectional surface, defines two neighborhoods of data elements around the two cusps 
in the initial cross-sectional surface and projects the neighborhoods onto the second cross- 
sectional surface. 

1 83. The program of claim 1 82, wherein each neighborhood projected onto the second 
cross-sectional surface includes data elements representing portions of the tooth and data 
elements representing the gum tissue surrounding the tooth. 
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1 84. The program of claim 1 83, wherein the data elements representing the tooth 
include voxels of one color and the data elements representing the gum tissue include voxels of 
another color. 

5 

1 85. The program of claim 1 83, wherein the computer, in locating the cusps in the 
second cross-sectional surface, locates the pair of data elements representing gum tissue that lie 
closest together, where each of the two neighborhoods projected onto the second cross-sectional 
surface includes one of the data elements in the pair. 

10 

1 86. The program of claim 1 70, wherein the computer, in applying the test to identify 
data elements on the gingival boundary, creates a series of roughly parallel 2D planes, each 
intersecting the dentition roughly perpendicular to an occlusal plane of the dentition, and each 
including data elements that represent a cross-sectional surface of the dentition. 

15 

1 87. The program of claim 1 86, wherein the cross-sectional surface in each 2D plane 
includes two cusps that roughly identify the locations of the gingival boundary. 

188. The program of claim 1 87, wherein the computer identifies the cusps in each 
20 cross-sectional surface. 

1 89. The program of claim 1 88, wherein the computer, in identifying the cusps, locates 
the cusps in one of the planes and then confines the search for cusps in an adjacent plane to a 
predetermined area in the vicinity of the identified cusps. 

25 

190. The program of claim 1 87, wherein the computer allows a human user to select 
data elements that roughly identify the locations of the cusps in a selected one of the cross- 
sectional areas. 
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191. The program of claim 188, wherein the computer searches for cusps in the 
selected cross-sectional area and confines the search to data elements lying within a 
predetermined area in the vicinity of the data elements selected by the human user. 

5 1 92. The program of claim 1 88, wherein the computer searches for cusps in an 

adjacent cross-sectional area and confines the search to data elements lying within a 
predetermined area in the vicinity of the data elements selected by the human user. 

193. A computer-implemented method for use in manipulating a digital model of a 
10 patient's dentition, the method comprising: 

obtaining a three-dimensional (3D) digital model of the patient's dentition; and 
analyzing the dentition model to determine the orientation of at least one axis of the 
dentition model automatically. 

15 1 94. The method of claim 1 93 , further comprising creating an Oriented Bounding Box 

(OBB) around the dentition model. 

195. The method of claim 194, wherein the dentition model has a z-axis that extends in 
a direction in which the OBB has minimum thickness. 

20 

1 96. The method of claim 1 95, wherein the z-axis extends from a bottom surface of the 
dentition model to a top surface of the model, and wherein the method includes automatically 
identifying the top and bottom surfaces of the dentition model. 

25 1 97. The method of claim 1 96, wherein one of the surfaces is substantially flat and 

another of the surfaces is textured, and wherein identifying the top and bottom surfaces includes: 
creating one or more planes that are roughly normal to the z-axis; and 
creating line segments that extend between the one or more planes and the top and 
bottom surfaces of the dentition model. 

30 
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1 98. The method of claim 197, wherein identifying the top and bottom surfaces 
includes identifying the surface for which all of the line segments are of one length as being the 
flat surface. 

5 1 99. The method of claim 1 97, wherein identifying the top and bottom surfaces 

includes identifying the surface for which the line segments have varying lengths as being the 
textured surface. 

200. The method of claim 193, wherein analyzing the dentition model includes 

1 0 selecting a two-dimensional (2D) plane that contains the axis and an arch-shaped cross section of 
the dentition model and identifying the orientation of the axis in this plane. 

201 . The method of claim 200, wherein the arch-shaped cross section is roughly 
symmetrical about the axis. 

15 

202. The method of claim 201, wherein analyzing the dentition model includes: 
identifying a point at each end of the arch-shaped cross section; 

creating a line segment that extends between the identified points; and 

identifying the orientation of the axis as being roughly perpendicular to the line segment 

20 

203. The method of claim 202, wherein identifying a point at each end of the arch 
includes: 

selecting a point that lies within an area surrounded by the arch-shaped cross section; 
creating a line segment that extends between the selected point and n edge of the 2D 

25 plane; 

sweeping the line segment in a circular manner around the selected point; and 
identifying points at the ends of the arch-shaped cross section at which the sweeping line 

segment begins intersecting the cross section of the dentition model and stops intersecting the 

cross section of the dentition model. 

30 
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204. The method of claim 201, wherein analyzing the dentition model includes 
identifying the orientation of another axis that is roughly perpendicular to the identified axis. 

205. A computer-implemented method for use in creating a digital model of an 
5 individual component of a patient's dentition, the method comprising: 

obtaining a 3D digital model of the patient's dentition; 

identifying points in the dentition model that lie on an inter-proximal margin between 
adjacent teeth in the patient's dentition; and 

using the identified points to create a cutting surface for use in separating portions of the 
10 dentition model representing the adjacent teeth. 

206. The method of claim 205, further comprising displaying 2D cross sections of the 
dentition model and receiving input from a human operator identifying approximate points at 
which the interproximal margin between the adjacent teeth meets gingival tissue. 

15 

207. The method of claim 206, wherein the dentition model includes a 3D volumetric 
model of the patient's dentition and the input from the human operator identifies two voxels in 
the volumetric model. 

20 208. The method of claim 207, further comprising defining a neighborhood of voxels 

around each of the two voxels identified by the human operator, where each neighborhood 
includes voxels representing the dentition model and voxels representing a background image. 

209. The method of claim 208, further comprising applying a computer-implemented 
25 test to select a pair of voxels, both representing the background image, that lie closest together, 

where each neighborhood contains one of the voxels. 

210. The method of claim 207, further comprising automatically identifying voxels on 
another 2D cross section that represent the interproximal margin. 

30 
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21 1 . The method of claim 210, wherein automatically identifying voxels on another 2D 
cross section includes: 

defining a neighborhood of voxels around each of the selected voxels, where each 
neighborhood includes voxels representing the dentition model and voxels representing a 
5 background image; 

projecting the neighborhoods onto the other 2D cross section; and 

selecting two voxels in the projected neighborhoods that represent the inter-proximal 
margin. 

10 212. The method of claim 211, wherein selecting two voxels in the projected 

neighborhoods includes selecting a pair of voxels, both representing the background image, that 
lie closest together, where each of the neighborhoods contains one of the voxels. 

213. A computer-implemented test for use in creating a digital model of an individual 
15 component of a patient's dentition, the method comprising: 

displaying an image of a dentition model; 

receiving input from a human operator identifying points in the image representing a 
gingival line at which a tooth in the dentition model meets gingival tissue; and 

using the identified points to create a cutting surface for use in separating the tooth from 
20 the gingival tissue in the dentition model. 

214. The method of claim 213, wherein the cutting surface extends roughly 
perpendicular to an occlusal plane in the dentition model. 

25 215. The method of claim 214, wherein creating the cutting surface includes projecting 

*at least a portion of the gingival line onto a plane that is roughly parallel to the occlusal plane. 

216. The method of claim 215, wherein creating the surface includes creating a surface 
that connects the gingival line to the projection. 

30 . ' 
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217. The method of claim 2 15, further comprising creating the plane by fitting the 
plane among the points on the gingival line. 

218. The method of claim 2 1 7, further comprising shifting the plane away from the 
tooth in a direction that is roughly normal to the plane. 

219. The method of claim 218, wherein shifting the plane includes creating a line 
segment that includes a point near the center of the tooth and that is roughly perpendicular to the 
plane. 

220. The method of claim 219, wherein the length of the line segment is approximately 
equal to the length of a tooth root 

22 1 . The method of claim 2 1 9, further comprising creating a sphere that has a radius 
equal to the length of the line segment and that is centered on the point near the center of the 
tooth. 

222. The method of claim 22 1 , wherein shifting the plane includes moving the plane 
along the line segment so that the plane is tangential to the sphere. 

223. The method of claim 222, further comprising receiving instructions from a human 
operator to slide the plane to a new position along the sphere. 

224. The method of claim 213, wherein the cutting surface extends roughly parallel to 
an occlusal plane in the dentition model. 

225. The method of claim 224, wherein the input received from the human operator 
identifies points that form two 3D curves representing gingival lines at which teeth in the 
dentition model meet gum tissue on both the buccal and lingual sides of the dentition model. 
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226. The method of claim 225, wherein creating the cutting surface includes fitting a 
surface among the points lying on the two curves. 

227. The method of claim 225, wherein creating the surface includes, for each tooth, 

5 identifying a point lying between the two curves and creating surface triangles having vertices at 
the identified point and at points on the two curves. 

228. The method of claim 227, wherein identifying the point includes averaging, for 
each tooth, x, y and z coordinate values of the points on portions of the two curves adjacent to 

10 the tooth. 

229. The method of claim 225, further comprising creating a surface that represents 
tooth roots. 

15 230. The method of claim 229, wherein creating the surface representing tooth roots 

includes projecting points onto a plane that is roughly parallel to the occlusal plane. 

23 1 . The method of claim 230, wherein creating the surface includes connecting points 
on the two curves to the projected points. 

20 

232. The method of claim 23 1, further comprising using the surface to separate 
portions of the dentition model representing the tooth roots from portions representing gingival 
tissue. 

25 233 . The method of claim 232, further comprising connecting the portions of the 

dentition model representing the tooth roots to the portion representing the tooth. 

234. A computer program, stored on a tangible storage medium, for use in 
manipulating a digital model of a patient's dentition, the program comprising executable 
30 instructions that, when executed by a computer, cause the computer to: 

obtain a three-dimensional (3D) digital model of the patient's dentition; and 
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analyze the dentition model to determine the orientation of at least one axis of the 
dentition model automatically. 

235. The program of claim 234, wherein the computer creates an Oriented Bounding 
5 Box (OBB) around the dentition model. 

236. The program of claim 235, wherein the dentition model has a z-axis that extends 
in a direction in which the OBB has minimum thickness. 

10 237. The program of claim 236, wherein the z-axis extends from a bottom surface of 

the dentition model to a top surface of the model, and wherein the computer automatically 
identifies^ top and bottom surfaces of the dentition model. 

238. The program of claim 237, wherein one of the surfaces is substantially flat and 

1 5 another of the surfaces is textured, and wherein, in identifying the top and bottom surfaces, the 
computer: 

creates one or more planes that are roughly normal to the z-axis; and 
creates line segments that extend between the one or more planes and the top and bottom 
surfaces of the dentition model. 

20 

239. The program of claim 238, wherein, in identifying the top and bottom surfaces, 
the computer identifies the surface for which all of the line segments are of one length as being 
the flat surface. 

25 240. The program of claim 238, wherein, in identifying the top and bottom surfaces, 

the computer identifies the surface for which the line segments have varying lengths as being the 
textured surface. 

241. The program of claim 234, wherein, in analyzing the dentition model, the 
30 computer selects a two-dimensional (2D) plane that contains the axis and an arch-shaped cross 
section of the dentition model and identifying the orientation of the axis in this plane. 
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242. The program of claim 241, wherein the arch-shaped cross section is roughly 
symmetrical about the axis. 

5 243 . The program of claim 242, wherein, in analyzing the dentition model, the 

computer 

identifies a point at each end of the arch-shaped cross section; 

creates a line segment that extends between the identified points; and 

identifies the orientation of the axis as being roughly perpendicular to the line segment. 

10 

244. The program of claim 243, wherein, in identifying a point at each end of the arch, 
the computer 

selects a point that lies within an area surrounded by the arch-shaped cross section; 
creates a line segment that extends between the selected point and n edge of the 2D plane; 
15 sweeps the line segment in a circular manner around the selected point; and 

identifies points at the ends of the arch-shaped cross section at which the sweeping line 
segment begins intersecting the cross section of the dentition model and stops intersecting the 
cross section of the dentition model. 

20 245. The program of claim 242, wherein, in analyzing the dentition model, the 

computer identifies the orientation of another axis that is roughly perpendicular to the identified 
axis. 

246. A computer program, stored on a tangible storage medium, for use in creating a 
25 digital model of an individual component of a patient's dentition, the program comprising 
executable instructions that, when executed by a computer, cause the computer to: 
obtain a 3D digital model of the patient's dentition; 

identify points in the dentition model that he on an inter-proximal margin between 
adjacent teeth in the patient's dentition; and 
30 use the identified points to create a cutting surface for use in separating portions of the 

dentition model representing the adjacent teeth. 
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247. The program of claim 246, wherein the computer displays 2D cross sections of 
the dentition model and receives input from a human operator identifying approximate points at 
which the interproximal margin between the adjacent teeth meets gingival tissue. 

5 

248. The program of claim 247, wherein the dentition model includes a 3D volumetric 
model of the patient's dentition and the input from the human operator identifies two voxels in 10 
volumetric model. 

1 0 249. The program of claim 248, wherein the computer defines a neighborhood of 

voxels around each of the two voxels identified by the human operator, where each 
neighborhood includes voxels representing the dentition model and voxels representing a 
background image. 

1 5 250. The program of claim 249, wherein the computer automatically selects a pair of 

voxels, both representing the background image, that lie closest together, where each 
neighborhood contains one of the voxels. 

25 1 . The program of claim 248, wherein the computer automatically identifies voxels 
20 on another 2D cross section that represent the interproximal margin. 

252. The program of claim 25 1 , wherein, in automatically identifying voxels on 
another 2D cross section, the computer: 

defines a neighborhood of voxels around each of the selected voxels, where each 
25 neighborhood includes voxels representing the dentition model and voxels representing a 

background image; 

projects the neighborhoods onto the other 2D cross section; and 

selects two voxels in the projected neighborhoods that represent the inter-proximal 30 

margin. 

30 



58 



WO 00/19935 



PCT/US99/23532 



253. The program of claim 252, wherein, in selecting two voxels in the projected 
neighborhoods, the computer selects a pair of voxels, both representing the background image, 
that lie closest together, where each of the neighborhoods contains one of the voxels. 

5 254. A computer program, stored on a tangible storage medium, for use in creating a 

digital model of an individual component of a patient's dentition, the program comprising 
executable instructions that, when executed by a computer, cause the computer to: 
display an image of a dentition model; 

receive input from a human operator identifying points in the image representing a . 
10 gingival line at which a tooth in the dentition model meets gingival tissue; and 

use the identified points to create a cutting surface for use in separating the tooth from the 
gingival tissue in the dentition model. 

255. The program of claim 254, wherein the cutting surface extends roughly 
15 perpendicular to an occlusal plane in the dentition model. 

256. The program of claim 255, wherein, in creating the cutting surface, the computer 
projects at least a portion of the gingival line onto a plane that is roughly parallel to the occlusal 
plane. 

20 

257. The program of claim 256, wherein, in creating the surface, the computer creates 
a surface that connects the gingival line to the projection. 

258. The program of claim 256, wherein the computer creates the plane by fitting the 
25 plane among the points on the gingival line. 

259. The program of claim 258, wherein the computer shifts the plane away from the 
tooth in a direction that is roughly normal to the plane. 
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260. The program of claim 259, wherein, in shifting the plane, the computer creates a 
line segment that includes a point near the center of the tooth and that is roughly perpendicular to 
the plane. 

5 26 1 . The program of claim 260, wherein the length of the line segment is 

approximately equal to the length of a tooth root 

262. The program of claim 260, wherein the computer creates a sphere that has a 
radius equal to the length of the line segment and that is centered on the point near the center of 

10 the tooth. 

263. The program of claim 262, wherein, in shifting the plane, the computer moves the 
plane along the line segment so that the plane is tangential to the sphere. 

15 264. The program of claim 263, wherein the computer receives instructions from a 

human operator to slide the plane to a new position along the sphere. 

265. The program of claim 264, wherein the cutting surface extends roughly parallel to 
an occlusal plane in the dentition model. 

20 

266. The program of claim 265, wherein the input received from the human operator 
identifies points that form two 3D curves representing gingival lines at which teeth in the 
dentition model meet gum tissue on both the buccal and lingual sides of the dentition model. 

25 267. The program of claim 266, wherein, in creating the cutting surface, the computer 

fits a surface among the points lying on the two curves. 

268. The program of claim 269, wherein, in creating the surface, the computer, for 
each tooth, identifies a point lying between the two curves and creates surface triangles having 
30 vertices at the identified point and at points on the two curves. 
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269. The program of claim 268, wherein, in identifying the point, the computer 
averages, for each tooth, x, y and z coordinate values of the points on portions of the two curves 
adjacent to the tooth 

270. The program of claim 266, wherein the computer creates a surface that represents 
tooth roots. 

271 . The program of claim 270, wherein, in creating the surface representing tooth 
roots, the computer projects points onto a plane that is roughly parallel to the occlusal plane. 

272. The program of claim 271 , wherein, in creating the surface, the computer 
connects points on the two curves to the projected points. 

273. The program of claim 272, wherein the computer uses the surface to separate 
portions of the dentition model representing the tooth roots from portions representing gingival 
tissue. 

274. The program of claim 273, wherein the computer connects the portions of the 
dentition model representing the tooth roots to the portion representing the tooth. 
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